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£\j ■ In this paper, we present Mie coefficients of double-layered sphere and consider the scattering 

problem, including the topics on field distribution, electromagnetic cross section, extinction spectra 
as well as some potential peculiar properties arising from the presence of left-handed media. The 
formulation presented here can be easily generalized to cases of multiple-layered spheres. 
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I. INTRODUCTION 



More recently, a kind of artificial composite media ( the so-called left-handed media ) having a frequency band where 
the effective permittivity ( e ) and the effective permeability ( /i ) are simultaneously negative attracts attention of 
many researchers in various fields such as materials science, condensed matter physics, optics and electromagnetism 
[1-5]. Veselago first considered this peculiar medium and showed that it possesses a negative index of refraction [6]. 
It follows from the Maxwell's equations that in this medium the Poynting vector and wave vector of electromagnetic 
wave would be antiparallel, i. e., the vector k, the electric field E and the magnetic field H form a left-handed 
system; thus Veselago referred to such materials as "left-handed" , and correspondingly, the ordinary medium in which 
k, E and H for m a right-handed system may be termed the "right-handed" medium. Other authors call this class 
of materials "negative-index media (NIM)", "double negative media (DNM)" [5] and Veselago's media. It is readily 
y—i , verified that in such media having both e and fi negative, there exist a number of peculiar electromagnetic properties, 
for instance, many dramatically different propagation characteristics stem from the sign change of the group velocity, 
including reversals of both the Doppler shift and the Cherenkov radiation, anomalous refraction, modified spontaneous 
' emission rates and even reversals of radiation pressure to radiation tension [2]. In experiments, this artificial negative 
electric permittivity media may be obtained by using the array of long metallic wires (ALMWs), which simulates the 
plasma behavior at microwave frequencies, and the artificial negative magnetic permeability media may be built up by 
using small resonant metallic particles, e. g., the split ring resonators (SRRs), with very high magnetic polarizability 
[3,7-9]. 

The extinction properties of a sphere ( single-layered ) with negative permittivity and permeability is investigated by 
Ruppin. Since recently Wang and Asher developed a novel method to fabricate nanocomposite SiC>2 spheres ( ~ 100 
nm ) containing homogeneously dispersed Ag quantum dots (2^5 nm ) [10], which may has potential applications 
to design and fabrication of photonic crystals, it is believed that the absorption and transmittance of double-layered 
sphere deserves consideration. In this paper, both the Mie coefficients and the electromagnetic field distributions in 
the two-layered sphere containing left-handed media irradiated by a planar electromagnetic wave are presented. 

1 IT. METHODS 

ft ■ 

It is well known that in the absence of electromagnetic sources, the characteristic vectors such as E, B, D, H and 
Hertz vector in the isotropic homogeneous media agree with the same differential equation 

VV ■ C - V x V x C + fc 2 C = 0. (2.1) 
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The three independent vector solutions to the above equation is [11] 

L = VV>, M = V x aip, N=|vxM (2.2) 

k 

with a being a constant vector, where the scalar function tp satisfies V 2 ?/> + k 2 ip = 0. It is verified that the vector 
solutions M, N and L possess the following mathematical properties 

M = Lxa=-VxN, L • M = 0, VxL = 0, V • L = V 2 ip = —k 2 ip, V • M = 0, V • N = 0. (2.3) 
k 

Set M = mexp(— iwt) and N = nexp(— iut), the vector wave functions M and N can be expressed in terms of the 
following spherical vector wave functions [11] 

m= m „ = T-^n(fcr)P„ m (cos0)->0i 2 - z n {kr)—^j- '-^m<t>h, 

sin u ou 

"F j^^KWI/T^IC^s. (2.4) 
In what follows we treat the Mie coefficients of double-layered sphere irradiated by a plane wave. 



A. Definitions 



We consider a double-layered sphere with interior radius cl\ and external radius a 2 having relative permittivity 
(permeability) ei (pi) and e 2 ( fJ-2), respectively, placed in a medium having the relative permittivity eo an d perme- 
ability fj,Q. Suppose that the double-layered sphere irradiated by the following plane wave with the electric amplitude 
Eo along the z-direction of Cartesian coordinate system [11] 

Ei = a x E CKp{ik Q z - iut) = E Q exp(-iwt) j " n fa + 1) ( m ii ) n ~ ir 4il) > 

71=1 

Hi = & y —E exp(ik z - iujt) = — — E exp(-iwt) V »" 7^-771 fm^ + m^) , (2.5) 



where 



mill - ±-^jn(kor)pi(cosO)Zcl>h - j n {k Q r)^™<PH, 
4i« = ! ^^j4fcor)P„ 1 (cos0)^^i 1 + -i- [fc onn (fc r)]' ^ >i 2 ± j^—^ [korj n (k r)}' P^(cos6l)» (2.6) 

with the primes denoting differentiation with respect to their arguments and fco being ^/eo^o^- In the region r < d\, 
the wave function is expanded as the following series 

00 2+1 

E t = E exp(-i W t) ]T *" ra( " + 1) (Onmoln " K^eln) , 
i 00 9 _i_ 1 

Ht = --^E exp(-ia;f) 2 i-^^y (^m^ + i<£n£ n ) . (2.7) 

^ n=l ^ ' 

Note that here the propagation constant in and n^ n which have been defined in (2.6) is replaced with k\, i.e., 

V e iMi^- I n the region a x < r < 02, one may expand the electromagnetic wave amplitude as 

00 Oil 

E m = E exp(-^i) £ (°n m iin - i^n^ + <m^ n - ib™n<$ n ) , 



.i(n+l) 

n=l v ' 

H m = --^2*> exp(-iart) ]T i"^Ty (Cm« + <n« + t%m$> n + <n^) . (2.8) 

^ n=l ^ ' 
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(3) (3) 

Note that here in order to obtain m y el ' n and n y ol ' n , the spherical Bcsscl functions j n in the spherical vector function 
m^ n and nj^ n is replaced by the Hankel functions , namely, the explicit expressions for and n^ n can also 



be obtained from (2.6), so long as we replace the spherical Bessel functions j n in (2.6) with the Hankel functions h 
Apparently, here the propagation constant in and n^ n should be replaced with k 2 , i.e., ^62^2^- 

In the region r > 0,2 the reflected wave is written [11] 



(i) 



^ 2tt, -|- l / 

E r = E exp(-iwi) V i n — — — I 
^ n(n + 1 V 



^ n(n + l) 

n=l ' 

H r = E cxp(-zwi) ) i n — — — 

Po^ ^ n(n + l) 



(3) 
a n m oln 



(3) , ■ r (3) 



(3) (3) 

where the propagation vector in m K el ' n and n y ol J n is replaced by k 0l i.e., ^/eo^o^- 



(2.9) 



B. Boundary conditions 

At the boundary r = a\, the boundary condition is given as follows 

ii x E m = ii x E t , ii x H m = ii x H t , 



where the unit vectors of Cartesian coordinate system agree with ii x 12 = 13, i 2 x 13 = ii, 13 x ii = 12- 
It follows from ii x E m = ii x E t that 

<%j n (N 2Pl ) + a™h^(N 2Pl ) = alj n {N lPl ), 
Nib™ [N2Pijn(N2Pi)}' + \N 2 p 1 hl 1 \N2p 1 )\ ' = N 2 b t n [N lPl j n (N lPl )}' , 



where p\ — fcoai, iVi — 



feo ' 



In the similar manner, it follows from ii x H m — \\ x H t that 

p 1 (a- {N2p 1 j n (N 2 p 1 )}' + < \N 2Pl h^ (N 2P i) \ = p 2 {a* [N 1 p 1 j n (N 1 p 1 )}'} 



N2P1 [b™j n (N2Pi) + b™h£\N2 Pl )] = N 1 p 2 b t Jn{N 1 pi)- 

At the boundary r = a 2 , in the same fashion, it follows from ii x (E; + E r ) = ii x E m that 

j n (p 2 ) + a I n h^(p 2 ) = a™j n (N 2 p2) + a™hW(N 2 p2), 
N 2 [P23n(p2)}' + N 2 V n [P2h£\p2)]' = C [N2P2MN2P2)}' + K [N2P2h^(N 2 p 2 ) 

where p 2 = fc a 2 , N 2 = ||. 

In the meanwhile, it follows from ii x (H; + H r ) = ii x H m that 

P2h^(p 2 )\ ' + lP2jn(P2)}'} = Pa f<C [N 2 p2j n (N2P2)}' + < \N 2 p 2 h^ (N2P2) 



p 2 i a n 



P2b T n h£\p 2 ) + P2UP2) = N 2 p a bZUN 2 p2) + N 2 pob™hi 1 \N2p 2 ). 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



C. Calculation of Mie coefficients 

1. Mie coefficient a r n 

According to the above boundary conditions, one can arrive at the following matrix equation 
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-jn(N lPl ) Jn(N 2Pl ) h£\N 2Pl ) 

-M2 [Nl P lj n (N lPl )]' Vl[N 2P ljn(N 2P l)}' Ml N 2Pl h { n ] (N 2Pl ) 

jn(^2P2) h£\N 2P2 ) 

Mo [N 2P2 j n (N 2P2 )}' \N 2P2 h { n 1] \N 2P2 ) 



"M2 






-h£\ P2 ) 

P 2hn\f>2) 



For convenience, the 4x4 matrix on the left handed side of (2.15) is denoted by 

/A n A 12 A 13 \ 

A 21 A 22 A 23 

A 32 A 33 A 34 

\ A 42 A i3 A 4i 



whose determinant is 



detA = A U A 22 A 33 A 44 - A U A 22 A 34 A 43 - AnA 32 A 23 A 44 + A 11 A 42 A 23 A 34: 
- A 21 A 12 A 33 A 44 + A 21 A 12 A 3i A 43 + A 21 A 32 A 13 A 44 - A 21 A 42 A 13 A 34 . 



So, the Mie coefficient a r n is of the form 



detA r 
dctA 



with 



A T 



/A n 


A 12 


A 13 


\ 


A 2 i 


A 22 


A 23 








A 32 


A 33 


3n{ P 2) 


V o 


A 42 


A 43 


M2 [ P 2jn( P 2)]' J 



whose determinant is 






V M2 [P 2 jn(p 2 )]' / 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



detA r = AnA 22 A 33 \i2 \p23n{Pi)\ - A n A 22 j n ( p2 )A 43 - AnA 32 A 23P2 [ P2 j n ( p2 )]' + AuA 42 A 23 j n ( P2 ) 
- A 21 A 12 A 33 fi2 [P2.7«(p 2 )]' + A 21 A 12 j n ( P2 )A 43 + A 21 A 32 A 13P2 [ P2 j n ( P2 )]' - A 21 A 42 A 13 j n ( P2 ). 



(2.20) 



2. Mie coefficient ¥ n 

According to the above boundary conditions, one can arrive at the following matrix equation 



( -N 2 [N lPl j n {N lPl )]' N, [N 2Pl j n {N 2Pl )\ N, [N 2Pl h^ (N 2Pl ) 
-N lP2 j n (N lPl ) N 2 ^j n (N 2Pl ) N 2Pjl h£ ] (N 2Pl ) 



V 



[N 2P2 j n (N 2P2 )}' 

N 2P0 jn(N 2P2 ) 



N 2P2 h£\N 2P2 ) 

N 2Pj0 h ( n\N 2P2 ) 



-No 






P 2h£\ P2 ) 



\ 







I 




K 






K 






\ b nJ 





For convenience, the 4x4 matrix on the left handed side of (2.21) is denoted by 



B 



/ Bn Bi 2 B\ 3 \ 

B 2 \ B 2 2 B 23 

B 3 2 B 33 B 34 

\ B 42 B 43 B 44 J 



\ 



N 2 [ P2 j n ( P 2)}' 
M2jn(p2) 

(2.21) 



(2.22) 



whose determinant is 
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deti? — B11B22B33B44 — -B11-B22.B34B43 — B11B32B23B44 + -B11-B42-B23-B34 
— -B21-B12-B33-B44 + -B21-B12-B34S43 + -B21-B32-B13-B44 — -B21-B42-B13S34. 



So, the Mie coefficient b T n is of the form 



bi = 



detS r 
dctB 



with 



/ B n 


B\2 


B13 







B21 


B22 


B23 










B32 


B33 


^2 [P2jn(P2)]' 




V 


B42 


B43 


P23n{p2) 


1 



B T 



whose determinant is 

detB r = BnB 2 2B33^2jn(P2) - BnB 2 2N 2 [p23n{p2)\ B43 - B n B 3 2B23P2jn(P2) 
- B2lB 12 B33fJ,2jn(P2) + B2lB 12 N 2 [p2jn(p2)}' B43 + B 2 lB3 2 B 1 3H23n{P2) 



B11B42B23N2 \p23n{P2)]' 
B2lB i2 Bi3N 2 [p2jn(P2)}' ■ 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



III. REDUCED TO THE CASE OF SINGLE-LAYERED SPHERE 



In order to sec whether the above Mie coefficients in the case of double-layered sphere is correct or not, we consider 
the reduction problem of double-layered case to the single-layered one when the following reduction conditions are 
satisfied: An = -A ±2 , A 2 \ = -A22 N± = N 2 , p\ = P2, Pi = P2- 

By lengthy calculation, we obtain 



detA = {A n A 2 3 - A 21 A 13 ) (A42A34 - ^32^44) 

= (A11A23 - A 2 lA 13 ) I jn{N 2 p2)p2 \p2h^ ) (p2)} - PO [N2P23n{N 2 p2^ (p 2 ) 



and 



Thus 



detA r = {AnA 2 3 - A 21 A X 3) {p [N 2 p23n(N2p2)]' 3n{p2) ~ in{N 2 P2)P2 [P2jn(P2)}'} 
r _. PO [N 2 p23n(N2P2)}' 3n(p2) - 3n{N 2 p2)P2 [p23n{p2^ 



3n{N 2 p2)P2 P2h i n\p2) ~ PO [N 2 p23n (N2P2)]' ' (P2 ) 



'!,(!)/ 



(3.1) 



(3.2) 



(3.3) 



which is just the Mie coefficient a r n of single-layered sphere [11]. 

In the same fashion, when the reduction conditions Bn = —B12, B 2 \ = — -B22 N\ = N 2 , p\ = p 2 , pi = P2 are 
satisfied, one can arrive at 



det-B = (-B11-B23 — B21B13) (B42B34 — B32B44) 

= (SllB23 - B2lBl3) l[N2P23n(N2 P 2)}' P2hl 1) ( P 2) ~ N^ j n {N 2P 2) L/^O^) 



and 



So, 



dctS r = (B11S23 - B2lBl3) {Nip 3n{N 2 p2) [P23n{p2)]' ~ [N 2 P23n (JV 2 p 2 )] ' p 2 j„ (p 2 ) } • 
Nip jn(N 2 p2) [p2jn(P2)}' - [A^jn (N2P2)] ' P23n{p2) 



bl = 



[N2p2jn(N 2 p2)}' P2h { n\p 2 ) ~ NlfX j n (N 2 p 2 ) [M^O^) 

which is just the Mie coefficient ¥ n of single-layered sphere [11]. 

This, therefore, means that the Mie coefficients of double-layered sphere presented here is right. 



(3.4) 



(3.5) 



(3.6) 



5 



IV. ALL THE MIE COEFFICIENTS INSIDE THE DOUBLE-LAYERED SPHERE 



A. Mie coefficients a* n , a„ and a™ 



The 4x4 matrix A t is 



( A 12 A 13 \ 

A 22 A 23 

3n{P2] A 32 A 33 A 34 

\ P-2 [P2jn(/f2)]' -442 A 43 A 44 / 



whose determinant is 

detA t = j n (P2)A 12 A 2 3A44 - j n (p 2 )A 2 2A 13 A44 - p 2 [p 2 j n {P2)}' A 12 A 23 A 34 + p 2 [P2Jti(P2)]' ^22^13^34- 

So, 

t detA t 



" dctA ' 



The 4x4 matrix A m is 



(An 





A 13 


° \ 


A 2l 





A 23 








3n{P2) 


A 33 


A 34 


V o 


M2 [P2j'n(P2)]' 


A 43 


A 44 J 



whose determinant is 

detA m = -A u j n (p 2 )A 23 A 44 + A lx p 2 [p 2 j n (p 2 )]' A 23 A 3i + A 21 j n (p 2 )A 13 A Ai - A 21 p 2 [p 2 j n (p 2 )}' A 13 A 3i . 

So, 

detA m 



detA ' 



The 4x4 matrix A™ is 



Am — 



/An Ai 2 \ 

A 2 i A 22 

A 32 jn{pl) A34 

\ A 42 p 2 [p 2 j n {p 2 )}' A44 / 



whose determinant is 

detAa = AuA 22 j n (p 2 )A 4 4 - AiiA 22 A 34 /i 2 [p 2 j n (p 2 )}' - A 2 iA 12 j n (p 2 )A 44: + A 21 A 12 A 3i p 2 [p 2 j n (p 2 )}' ■ 

So, 

m _ dctAg! 
° n " "detA"' 



It is readily verified that under the reduction conditions An = —A\ 2 , A 2 \ = — A 22 , N\ = N 2 , p\ = p 2 . 
one can arrive at 

detAt = detA m , = a™, dctA ffl = 0, a* = 0, 
which means that the above Mie coefficient can be reduced to the those of single-layered sphere. 



B. Mie coefficients 6™ and 6™ 



The 4x4 matrix B f is 



B t = 



( Bia Bis \ 

-B22 B23 

^2 [p2jn(P2)]' #32 B33 -834 

V Vl3n(pl) B A2 -B43 B44 / 



(4.11) 



whose determinant is 

detBt = ^2 [p2jn(P2)]' B 12 B 23 B 44 - ^2 [p2j'n(P2)]' B 22 B 13 B 44 - p 2 in{P2)B 12 B 23 B 34 + p 2 jn(P2)B 22 B 13 B 34 . (4.12) 

So, 

, t detB t 



dctB 



(4.13) 



The 4x4 matrix _B m is 



B 2 i 



B13 \ 

B 23 

N 2 [p 2 j n (p 2 )}' B 33 B 3i 

V P2jn(P2) -B43 #44 / 



(4.14) 



whose determinant is 

detB m = -B U N 2 [p 2 jn(p2)}' B 23 B 44 + Bnp 2 j n (p 2 )B 23 B 34 + B 21 N 2 [p 2 jn(p 2 )]' B13S44 - B 21 p 2 j n (p 2 )B 13 B 34 . 

(4.15) 

So, 

detB m 



dctB 



(4.16) 



The 4x4 matrix B r „ is 



Brn = 



/ B u B 12 \ 

£>2i S22 

B 32 N 2 [p 2jn {p 2 )]' B 34 

V B 42 A*2in(P2) s 44 y 



(4.17) 



whose determinant is 

detB ffl = BuB 22 N 2 [p 2 j n (p 2 )}' B 44 - B Y1 B 22 B 34 p 2 j n {p 2 ) - B 21 B 12 N 2 [p 2 j n (p 2 )}' B 44 + B 21 B 12 B 34 p 2 j n (p 2 ). (4.18) 
So, 

det-Bm 



b m = 



dctB 



(4.19) 



It is readily verified that under the reduction conditions Bu = — B\ 2 , B 2 \ = —B 22 , N± — N 2l p\ = p 2 , pi = p 2 , 
one can arrive at 

detB t = detB m , b\ = 1%, dct£? ffl = 0, b™ = 0, (4.20) 
which means that the above Mie coefficient can be reduced to the those of single-layered sphere. 
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V. DISCUSSION OF SOME TYPICAL CASES WITH LEFT-HANDED MEDIA INVOLVED 



In this section we briefly discuss several cases with left-handed media involved. 



(i) If the following conditions h„\N 2 pi) = 0, 
h ( n 1 \ P2 ) 



h { n\N 2Pl ) 



0, h { n 1) (N 2P2 ) 



h ( n\ P2 ), h> n L \N 2P2 ) 



N 2 = fi 2 = e 2 = —1 and p — +1 arc satisfied, then it is easily verified that 
det^4 = 0, dct_B = and ai = oo, bi 



CO. 



(5.1) 



Thus in this case the extinction cross section of the two-layered sphere containing left-handed media is rather large. 



(ii) If the following conditions h ( n ) (N 2 pi) = 0, [/^(JVapi)] = 0, h£ ] \N 2 p 2 ) = j n (p 2 ), [h { n\N 2 p 2 ) 
N 2 = p 2 = e 2 = —1 and ^ = +1 arc satisfied, then it is easily verified that 

detA r = 0, detB r = and a\ = 0, b T n = 0. 



[jn(P2 



(5.2) 



Thus in this case the extinction cross section of the two-layered sphere containing left-handed media is negligibly 
small. 

(iii) The case with iVi = pi = e\ = —1, N 2 = p 2 = e 2 = 0, N a = /i = eo = +1 is of physical interest, which 
deserves consideration by using the Mie coefficients presented above. 

Based on the calculation of Mie coefficients, one can treat the scattering problem of double- layered sphere irradiated 
by a plane wave. The scattering and absorption properties of double-layered sphere containing left-handed media can 
thus be discussed in detail, which are now under consideration and will be submitted elsewhere for the publication. 
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